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Abstract 

Using the Bargmann-Husimi representation of quantum mechanics on a torus 
phase space, we study analytically eigenstates of quantized cat maps ||9|. The linear¬ 
ity of these maps implies a close relationship between classically invariant sublattices 
on the one hand, and the patterns (or ‘constellations’) of Husimi zeros of certain 
quantum eigenstates on the other hand. For these states, the zero patterns are crys¬ 
tals on the torus. As a consequence, we can compute explicit families of eigenstates 
for which the zero patterns become uniformly distributed on the torus phase space 
in the limit h —> 0. This result constitutes a first rigorous example of semi-classical 
equidistribution for Husimi zeros of eigenstates in quantized one-dimensional chaotic 
systems |Q|. 


1 Introduction 

A major question in “quantum chaos”, i.e. the study of quantum systems for which 
the classical limit is chaotic, lies in the structure of eigenstates. Indeed, a ‘chaotic 
counterpart’ to the WKB-Ansatz, which holds in the case of an integrable system, 
is still missing (this Ansatz gives simple asymptotic formulas for individual eigen¬ 
states). The only proven results so far, collectively referred to as ‘Schnirelman’s 
theorem’, deal with the phase-space measures associated to eigenstates of a classi¬ 
cally ergodic system, in the semi-classical limit. 

That is, to any quantum state b?i(Q') there corresponds a phase-space measure 
density defined using coherent states, called the Husimi density of V’/i [@1- 

Now, if one quantizes an ergodic Hamiltonian system and considers a sequence of 
eigenstates {b?in}neN s.t. hn —> 0, with eigenvalues En —> E, then the associated 
Husimi measures will almost always converge to the Liouville measure pE over 
the energy surface that is, for every smooth observable f{q,p), the averages 
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If{q,p)H^M:P)dqdp tend to JY,^f{q,p) dpE [|, IJ]. When the ergodic 

classical system is not a flow, but rather a map, the same property holds for almost 
any semi-classical sequence of eigenstates of the quantized map, the Husimi density 
now converging towards the Liouville measure over the whole phase space [^, |^ . 

These results are very appealing in that they do not depend on the details 
of the system considered, apart from its classical ergodicity. On the other hand, 
Schnirelman’s theorem allows a minority of non-ergodic eigenstates (a set of asymp¬ 
totic density zero), and we ignore what these might look like: some of them could 
for instance show a strong concentration along an unstable periodic orbit (a ‘scar’), 
even in the semi-classical limit. Furthermore, the information provided by this the¬ 
orem remains measure-theoretical, which is much weaker than a WKB-type Ansatz. 
For this reason, we intend to study as precisely as possible the structure of eigen¬ 
states for a class of analytically tractable chaotic systems, the automorphisms of the 
2-torus, or generalizations of Arnold’s cat map |^]. The quantization of these highly 
ergodic maps was first performed by Hannay and Berry |p, and has been rederived 
in different ways since then @11, 


P3 ]; their spectral properties were analyzed 
by semi-classical means in [^, ^|, and eigenstates were studied analytically by 

11, sp¬ 
in this paper we exhibit a certain class of quantum cat eigenstates, directly 
within the Bargmann-Husimi representation of quantum mechanics on the torus. 
This phase-space formalism, together with the linearity of the classical maps, allow 
us to associate a classically invariant sublattice of the unit torus with a family of 
eigenstates of the quantum map. Precisely, we show that the Husimi functions 
of these states (for which we have analytical expressions) are periodic w.r. to the 
invariant sublattice. We call these states crystalline, due to the patterns formed by 
the zeros of Moreover, we can show that in the limit of a large number of sites, 
these invariant sublattices (and the crystals they support) become equidistributed 
over the torus. Our derivation is a first step towards a still missing proof of the 
‘uniform distribution of the Husimi zeros’, conjectured by Leboeuf and Voros |Q, ^ 
for one-dimensional chaotic eigenstates. This conjecture represents a more precise 
statement than Schnirelman’s theorem, inasmuch as the zeros completely encode 
the eigenstate. Unfortunately, they behave still unpredictably in the semi-classical 
limit, and are in general not easy to catch analytically. What saves us here is 
the linearity of the classical map, which yields explicit (and tractable) formulas for 
eigenstates. 

Before proceeding further, let us sketch the plan of the paper. After recalling 
the basics of quantum mechanics on the torus and describing the Bargmann and 
Husimi representations in this framework (section |^, we describe the family of 
dynamical systems to be studied, i.e. the quantized hyperbolic automorphisms of 
the 2-torus (section |3|) . After explaining how to build naturally a quantum eigenstate 
from a classical invariant sublattice (section |^, we try to classify these lattices 
as thoroughly as possible, following the seminal work of Percival and Vivaldi [“M] 
(section ^). We estimate the uniformity of these invariant sets in the semi-classical 
limit (section P) , and then explicitly build different types of eigenstates (sections 
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HI) associated to these sets; the uniformity results can then be applied to the zero 
patterns of their Husimi representations. 


2 Quantum mechanics on the torus 

The classical maps we study here are defined on the unit torus phase space T^. 
Before quantizing such a map, we shall first define the quantum kinematics, i.e. the 
Hilbert space(s) of quantum states corresponding to this classical phase space. We 
will not resort to geometric quantization, as was done in but use instead a more 
pedestrian approach, which yields basically the same results (it amounts to a choice 
of polarized sections parallel to the p axis). Namely, we build our quantum states by 
imposing the following quasi-periodicity conditions upon a one-dimensional wave- 
function i.e. a distribution in 5'(R) |Q, |2^ : 

{q + l\^P) = e2--i(g|V.) (1) 

{p+l\^|J) = {p\iP) 


with usual Dirac notations for position and momentum representations. The Flo- 
quet parameter k = ki +iK 2 can take any value on the dual torus. These conditions 
force Planck’s constant to take discrete values h = (27rA^)“^, with N a positive 
integer. For given N and k, the quasi-periodic states form a Wdimensional Hilbert 
space of distributions TCn,k, having for example the following basis: 

kj)iV,K = \qj +m), for j = 0, • • •, - 1 (2) 

m£Z 


where the allowed position peaks are discrete, qj = 27 rh{K 2 + j)- The Hermitian 
structure on 7iN,K is simply given by {qk\qj)N,K = for j, - 1. 

Through a finite Fourier transform of \qj)N,K, one obtains the basis built by pe¬ 
riodization of momentum eigenstates, \pm)N,K, with the dual discretization pm = 
2TTfl{Ki + m). 

It can be helpful to work in a phase-space representation, since we are interested 
in semi-classical features of states in TCn,k (which arise for large values of N). We 
will use the Bargmann representation of quantum states on the real line |j^], simply 
extending it to the above distribution spaces by linearity |Q. This transformation 
maps a Schrodinger wave-function {q\ip) into an entire function {z\^jJ), through the 
integral kernel 


{z\q) 






zq 


(3) 


where jz) is a (Weyl) coherent-state, which is localized at the phase-space point 
{Qo,Po) given hy z = ■ Although \z) depends on fi, this dependence will not 

be explicit in our subsequent notations : the bra {z\ will always correspond to the 
value oifi {= l/27rA') labelling the attached ket \' 4 ’)n,k- 
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The combination of the equations (||,^) leads to the following representation of 
the basis \qj)N,K '■ 

«3 (""(«- ‘ 9 - “ (4) 

where we used Whittaker and Watson’s [^] definition for the Jacobi theta function, 
9^{Z\t) = Y, for Z,t GC, 9(r) > 0. (5) 

nEZ 


It can be convenient to use a more general definition of theta function |W|: 


0Ar,g. (r, Z,t) = e 


— 2i7rW£ ^ ^ ^inNr'y'^—2i7TN'yZ 
'y^Qj +Z 

= g-2i^Aft g-2i^iVZg,. ^i^Nrq] (^Z - Tqj)\NT), 


( 6 ) 


in which the basis vectors of TCn,k read as 

- 9-’ ®9 “ y) ■ 

We will use the latter type of theta functions when studying the action of the 
symplectic group upon biperiodic states, in the next section. On the other hand, 
the analytic properties of the “old” theta function 6*3 are thoroughly described in 
p^ , and we need them too. 

In the Bargmann representation, the vectors of 7i are completely character¬ 
ized by the following quasi-periodicity properties : 

(^+l/^/2|V^)^^,. = e2--ie-^(V2+v^.) (8) 

{z + i/V2\il;)N,. = e2--2e-iV(i/2-iV2z) 


The functions {z\^jJ)]\^,K are thus quasi-periodic w.r. to the torus Tc of periods 
Xj they actually represent particular holomorphic sections of a complex 
line bundle over this torus \U, |3^ . 

The scalar product of two vectors in 7iN,K is easily computed in Bargmann’s 
representation, thanks to a closure formula: 


(V’IV’%,K = 2 [ dRzdX^z {z\'ilj)jq^{z\'il}')N,K e 

JTc 


(9) 


Meanwhile, the positions of the zeros of a function {z\'ip)]s[,K are well-defined on 
Tc; there are exactly N of them, and their sum is constrained by 


N /I A 

\/2 '^^zk = N ^j - mod [1,i]. 


( 10 ) 
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The knowledge of the zeros allows to rebuild the quantum state, up to a constant 
factor C, from the unique function xiz) = {z\q = 0) jv=i,k=o ^ if we know the zeros Zk 
of a Bargmann wavefunction {z\'ij:)N,K (which satisfy (|l^)), we recover this function 
multiplicatively: 

{z\'tlj)N,K = C exphiTz'^izk -+ ■ ( 11 ) 

(The whole product is independent of the determination of the zeros we have used 
to build it.) This factorization property makes the Bargmann zeros a particularly 
interesting tool to construct or characterize a quantum state |^, Q, especially in the 
semi-classical framework. On the one hand, the set of zeros forms an exact and 
minimal representation of the quantum state, on the other hand the zeros live in 
the classical phase space, as opposed to the usual Schrodinger coefficients. This is 
the reason why we intend to investigate their positions for eigenstates of quantized 
cat maps. This task is a priori not trivial because the quantum state depends non- 
linearly of its zeros (the difficulty is similar to that of relating the coefficients and 
the roots of a polynomial). 

We are however able to give a complete description of the Bargmann zeros for 
particular states, namely the basis states of 'Hn,k (0), since we know where 9^{ttZ\t) 
vanishes ||^ : it has a unique zero in the fundamental torus of modulus r, lying at 
Z = {\ + r)/2. Therefore, the N zeros of {z\qj)M,K have the positions 

V2zk = qj + ]^-i '^mod [l,i], fc = 0, (12) 

= ^(i^) + ^+ i^ mod [l,i], 

i.e., they lie along a vertical line on the torus, at maximal distance from the line 
{q = qj} (see figure 1). Actually, the Bargmann function {z\qj)]\[,K verifies the ‘fine’ 
quasi-periodicity property 

{z + i/NV2\q,)N,. = e-(V27V-iv/2.) , ( 13 ) 

To close this section, let us describe the phase-space measure we can build from 
the Bargmann representation, i.e. the Husimi measure density |Q]. For any state 
\^p)n in 5'(]R), it is defined by the positive density 

H^{z,z) = \{z\'ip)\‘^ e-f . (14) 

For a state \'tp)N,K, the Husimi density function is single-valued on Tc- Like the 
Bargmann function, it factorizes easily in terms of its zeros z^ and of the Husimi 
density of x(z): 

N 

z,z) = K Y[ H^{z +Zk)- (15) 

k=l ' 

where K = \C\^ is the constant which makes H^{z,z) the density of a probability 
measure on Tc, as long as \'iP)n,k is normalized (cf. eq. ^). 
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What is known about the Husimi density ? The factorization proper¬ 
ties of Bargmann and Husimi functions mean that all the information about them 
(except an arbitrary factor) is contained in the positions of their zeros. We now 
roughly sketch the relation between the global features of the Husimi densities on 
the one hand, the patterns of zeros on the other hand, for two different types of 
eigenstates. 

For an integrable system (e.g. a time-independent Hamiltonian system on the 
torus), the Husimi measures of a sequence of eigenstates {{z\'ti))N,K}NeN with energies 
Ei\i ^ E are known to concentrate along the orbit of energy H [^|. In the same 
time, the zeros of these Husimi densities line up along certain classically defined 
lines, far away from the classical orbit (namely, anti-Stokes lines : see for instance 
0). These results are mere consequences of the WKB asymptotic form of the 
eigenstates of integrable systems. As an example, one can consider the position 
eigenstates \qj)N,K (see figure 1); for any Q G [0,1], we can select a sequence of 
states {\qj(N))N,K}N£n, such that qj(N) —> Q as —> oo. The zeros of {z\qj(^N))N,K 
lie on vertical lines which converge to the line {q = Q + 1/2}, whereas the Husimi 
densities concentrate semi-classically upon the line {q = Q}- 


On the opposite, in the case of classically ergodic dynamics (e.g. some kicked 
systems, or the quantum cat maps we will study here), the only analytical re¬ 
sult concerning the eigenstates is provided by Schnirelman’s theorem, described in 
the introduction. For almost every sequence of eigenstates {V'TvjAreN of the corre¬ 
sponding quantum maps, the Husimi measures converge to the Liouville measure 
on in the weak-* sense, which means that for any continuous function / on the 
torus, fj 2 f{q-,p)H^^{q,p) dqdp converges to the ergodic average fj 2 f{q,p) dqdp as 
N ^ oo 1^, (the interested reader can find precisions about the weak-* topology 
in 0). 

This weak-* convergence does not provide very precise information about the 
eigenstates. Moreover, it would be convenient to obtain a direct description of the 
phase-space minimal data formed by the zeros. Numerical calculations have yielded 
the following conjecture ii: one expects the zeros of a chaotic eigenstate to spread 
uniformly over the whole phase space (i.e. the torus) further studies showed 

that the zeros locally behave as roots of random polynomials ||5| (in particular, they 
tend to repel each other at short distance). This uniformity would be, in a certain 
sense, dual to the uniformity of the Husimi measure itself. In the following, we 
will be able to prove such a semi-classical equidistribution for very particular cat 
eigenstates, in which case the zero patterns are crystalline. 


3 Classical and quantum cat maps 


The classical transformations we consider are the hyperbolic automorphisms of the 


unit torus phase space , given by matrices S = 


a b 
c d 


in SL(2, Z) with tr(S') = 


a -|- d > 2 (this has the useful consequence 6c / 0). Such a matrix acts linearly on 
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the torus; 


X = ^ ^ J mod 1 —!■ S X mod 1. (16) 

These diffeomorphisms of are known to be fully ergodic and mixing ; in partic¬ 
ular, for a sufficiently smooth observable on the torus f{x), we have the property: 


for almost all x G T^, 


1 ^ r 

n=l 


(17) 


where (Px is the usual Liouville measure on T^. The condition “for a. a 
gratuitous: in the present case, all points on the torus with rational coordinates are 
periodic under the action of S, so the ergodic property (p^ ) does not hold for them. 
The periodic orbit structure of classical cat maps was largely studied in I'M, ‘M, M\, 
and we will use some of those results in the following. 

The quantization of S can be performed in two steps First, any transfor¬ 
mation S in SL(2,IR) acting on the whole plane can be considered as the time-one 
map of a certain quadratic Hamiltonian flow, as long as tr(S') > —2. For arbitrary 
h, evolving the Weyl-quantized Hamiltonian during a unit time yields the unitary 
operator Us on L 2 (K), given by the integral kernel 


W\Us\q) = 


1 


\/2mfib 


eft 


Uidq'^-2q'q+ag^)/2b] 


(18) 


(This corresponds exactly to the Van Vleck quantization scheme, since the term 
between brackets is the generating function of the symplectic transformation S). 
The ambiguity of the square-root sign is natural: the quantization S —> Us only 
provides a representation of SL(2,R) up to a sign, called the metaplectic represen¬ 
tation Mp(2,K). Elements of the group Mp(2,]R) are precisely specified by both the 
symplectic matrix S they are built from, and a holomorphic function on the upper 
half-plane ji^r) s.t. = br + d. This function appears naturally if we write the 

above kernel in a mixed Bargmann-Schrodinger representation: 

(19) 

This formula is valid for any value of h. However, in the following, working on the 
torus imposes the values h = (27rV)“^. 

Now that we have quantized the symplectic transformation S G SL(2,R) into a 
unitary operator Us on T 2 (M), we extend the action of Us to the distribution spaces 
TCn,k by linearity (see Appendix A). Using the notation S{t) = we obtain 


{z\Us\qk)N,K 


{z\Us\qk) O3 (^TrN[S{i)qk 


-iV2' 

\\f2z 
lb + d 


z 

\b + d 


V’ 


+ -] 


S{i)N 
2 


Ki Z 
Qk - 1 — 

^ N 2 


( 20 ) 

6^2 \ 
ib + d) ■ 
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This represents a generalization (to all k) of the action of Mp(2,Z) on theta func¬ 
tions, described for instance in More precisely, the transformation Us acts onto 
the function 0Af,gs,(T Z — ki /N,t) according to: 






( 21 ) 


To obtain (^), i.e. the transformation rule for the basis states of 'Hn,k-, we only 
need to select the three parameters: r = i, Z = iV2z, t = -iz‘^/2 + qkKi/N. 

For an arbitrary transformation S G SL(2,R), the transformed Bargmann func¬ 
tion is not an element of it is a theta function, quasi-periodic w.r. to another 

lattice: 


{z + ^^\Us\i^)N,. = (z|C/5|V’)iV,. 

(z+^^^|C/5|V^)7V,k = e^Ar[l!^+V2.(-6-id)] 


( 22 ) 


If S is a matrix with integer entries, this lattice coincides with the period lattice of 
Tc by modularity. The space 7iN,K is then mapped by Us into 'Hn,k’-, with |^, 



= S 



mod 1. 


(23) 


For simplicity, we will generally use the same notation Us for the unitary operator 
on L 2 (M) and the unitary N x N matrices connecting TCn,k and 'Hn,k'j since the 
latter are extensions of the former to some distribution spaces. In order to avoid 
the inhomogeneous term in eq. (|^), we choose (cf.|P, ^]) to restrict ourselves to 
transformations S of the following types: 


odd even 
even odd 


or 


even odd 
odd even 


(24) 


These transformations form a subgroup of index 3 in SL(2,Z), named Fi ^2 in the 
literature [^]. For all values of N, the spaces 7iN,K for k = 0 (periodic wave- 
functions) and K = 0 = (1 -|- i)/2 (anti-periodic wavefunctions) are both invariant 
under all operators Us, S G ri, 2 - We thus obtain two families of (non-faithful) 
finite-dimensional representations of this group. On the other hand, according to 
equation (^) an individual operator Us can act as an endomorphism on 7iN,K for 
some other values of k as well (see Appendix A, eq. ®)- 


4 Strategy: from zeros to eigenstates 

From now on, we give ourselves a fixed hyperbolic transformation S G Fi^ 2 ; and 
study the eigenstates of the corresponding matrices Us on the above-mentioned 
invariant Hilbert spaces, for any value N. 












Analytical expressions have already been obtained for eigenstates of Us on TC]\f,o 


fig , 1^ , ^|. The following trick was used: for all N, there exists an integer period 
p{N), s.t. 1, where 1 is the identity over 7iN,o, and cr{N) is a 

phase. Then, starting from any state \'iP)n,o, suitable combinations of its images 
Up'ip)N,o for A; = 1, • • • ,p{N) yield an orthogonal family of eigenvectors spanning 
a subspace of with eigenangles regularly spaced on the unit circle. One can 

then repeat this construction from a new vector \'ip')N,o orthogonal to this subspace 
and so on, until a complete basis is formed. Actually, we will use a similar approach 
in section (^. 

Unfortunately, the resulting formulas for the eigenstates are not very sugges¬ 
tive of their ’’shape”, in particular the limit —> oo looks quite untractable. By 
contrast, the main task of this article is to construct subfamilies of eigenstates char¬ 
acterized by interesting phase-space features, namely crystalline structures, when 
expressed in the Bargmann representation. Our approach will be similar to the 


one developed by Degli Esposti, Graffi and Isola |^]: we will assign to N certain 
values, selected according to their arithmetical properties (we generalize the prime 
values considered in [pl| ), and for these N we will exhibit some eigenstates showing 
lattice-like, or more generally crystalline zero patterns. We will use the close link 
between quantum and classical motions, due to the linearity of the classical map. 
Indeed, our main tool is the following statement : 

Proposition 1 The action of Us on any basis state \qj)N,K transforms the zeros 
of its Bargmann function {z\qj)M,K according to the classical map S (up to complex 
conjugation). 


This statement follows directly from eq. (^), since we know where the function 03 
vanishes . Calling the transformed zeros let us relate the positions of Zk (eq. 12) 
and zi 


k ‘ 


V2 Zk = {qj+ 1/2) +i— {1/2 + k - Ki) mod[l,i] 

V2 z'k = {a - ic) {qj + 1/2) + {-b + id)—{1/2 + k - Ki) mod[l,i]. 
This can be written, with obvious notations : 


(25) 


Zk = 



Zk mod [l/\/2, i/\/2] 


-/ 

Zk 


= S Zk mod [l/\/2, i/\/2] 


(26) 


(the complex conjugation is due to the convention z = p^)- 

The above lemma can be extended to the images of the basis states \qj)N,K 
under Mp(2, Z), thanks to group properties. Indeed, if we note Zk,z'f^, z'/ the zeros of 
{z\qj)N,K, {z\Us'\qj)N,Ki {z\Us"Us'\qj)N,K respectively (indexing the zeros correctly), 
the group law Us>'Us' = +Us”s' combined with eq.(26) imply 


4 = 3” S' Zk = S''zl 


(27) 
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move classically under Us"- 

is only valid for very special states, 


which shows that the zeros of {z\Us'\qj) n, i 

We stress that the transformation law 
and is not at all generic in 'Hn,k (In section (7^), we build another type of states for 
which (p6[ ) also applies, but this construction is possible only when is a perfect 
sqnare). 

In the following, any state {ip) obtained as the image of a basis state \qj)N,K under 
an operator Us will be called a “pure theta state”, since its Bargmann function can 
be written as a single theta function times an exponential factor (eq. |^ ). According 
to eq.(^), the N zeros of such a state make up a discrete ‘string’ on the torus Tc, 
of the form 


Zk = Zg + k{n + im)/NV2, for A; = 0, • • •, — 1. 


(28) 


Conversely, a state for which the zero pattern takes the form (28), with coprime 
integers n,m, is a pure theta state in the space 7iN,K (/^ is drawn from the sum rule 
(H))- For each N, we thus have a characterization of the set of pure theta states 
{Us\qj)N,K, S G SL(2,Z), K G [0, Ip, j = 0,..., A" — 1} through their Bargmann 
zeros. 

Performing a dilation of factor N\P2 and a reflection w.r. to the real axis, we 
map the lattice of zeros into an affine lattice A = Vo + Ajin, where Ajin is the integer 

n 


linear lattice generated hy v = 


—m 


modulo the square torus of side N. From 


the above properties, we deduce the following key property: 


Proposition 2 If there exists a pure theta state \'tp)N,K whose corresponding affine 
string A is invariant under S, then this state is an eigenstate of Us (and automat¬ 
ically belongs to one of the spaces Un,k on which Us acts as an endomorphism). 
Besides, the invariance of A implies that of Aijn by linearity. 


In mathematics (see for instance |^|, p.81-84), the two-dimensional lattice of 
integer points modulo N constitutes the module (here Zn denotes the ring of 
integers modnlo N; in this context, the integer N is called a modulus |^^). The 
snblattice Aun is a submodule of type 1 (because it is generated by v alone in Z^v), or 
principal submodule, which is free if n, m, N are globally coprime, or, equivalently, 
if Alin contains N points (by definition, a module is free if its generators are linearly 
independant). The submodules of Z^ we will consider will be principal, unless 
stated otherwise. 

In view of the above proposition, our first task will thus consist in identifying the 
integer strings Ann invariant modulo N, which actually amounts to diagonalizing 
S over Z^. The next section deals with the classification of these free principal 
eigenmodnles, for a given transformation S G Fi ^2 and a general value N. In a 
following stage, we will go back to quantum mechanics, giving explicit formulas 
for the eigenstates corresponding to the invariant strings A translated from Aijn. 
Finally, to obtain crystalline eigenstates, we will use sublattices also built through 
translations of Aun, but not invariant under S. 
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5 Free eigenmodules 


Given a fixed transformation S, we classify its free eigenmodules A^n for a general 
value of N (in this section, we will drop the subscript un and denote a linear string 
A). This classification deeply relies on the arithmetical properties of N. The values 
N for which there exists at least one free eigenmodule are called admissible (this 
property obviously depends on the matrix S). We use here a different formalism 


from Percival and Vivaldi’s |24], but our results highly overlap with theirs (we do 
not need the complete classification of periodic orbits these authors have worked 
out). 

When diagonalizing S' as a hyperbolic matrix acting on we obtain the fol¬ 
lowing formulas for eigen-values and -vectors : 


A± = 


. + d± \/{a + dy — 4 . 


is associated to the eigenvector 


k± = 



(29) 


The important quantity here is the argument under the square-root, i.e. the dis¬ 
criminant of S. We choose to normalize this discriminant as D = (^^)^ — 1 in the 
following {D is integer according to eq.(|^)). For the case of hyperbolic maps we 
consider (|a -|- d| > 2 ), the square-root of D in R is always irrational, so there is no 
eigenvector of S on R^ with rational coordinates. 

On the other hand, S can have integer eigenvectors when acting on provided 
the discriminant D admits an integer square root modulo N. This property is an 
arithmetical one, and we will need to factorize N as a product of primary numbers 
(i.e. powers of a prime). 

Before proceeding further, we fix our arithmetical notations: if M is a divisor 
of N, we note M\N', the greatest common divisor (gcd) of a family of integers will 

be noted (ni n n 2 n, • • • ,r\nr), and Ni N/{N n 1) for any I G Zjv. The simplest 
case to deal with is A^ = P a prime number; Zp is then a field, and Zp a vector 
space 

the primes P > 3 


The existence problem of a square-root of D modulo P yields a partition of 
12|: if there exists an integer x coprime to P s.t. = D mod 


P, D is said to be a quadratic residue of P, or equivalently, P is said to be splitting 
w.r.to D (this implies in particular (P n P) = 1). The Legendre symbol 0,111 is 
then defined as: 


= 1 if P splits w.r. to D 


I = 0 if P is a divisor of D 


(30) 


D 
P 
D 
P 
D 
P 

(In the following, we might forget the “w.r.to P” when mentioning a splitting prime). 

The case where P is inert will not be considered since it yields no eigenmod¬ 
ule. We will treat the case of splitting primes in the main text, and defer two 
complementary cases (P|P and P = 2) to Appendix B. 


— I = —1 otherwise (P is said to be inert w.r.to P) 
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To go from prime to primary values N = P^, we will use the following lemma 

0 : 

Lemma 1 To any solution kp of an algebraic equation f{k) = 0 modulo P s.t. 
f'{kp) ^ 0 mod P, there corresponds for any (3 > 1 a unique solution k^ modulo 
N = Pl^; moreover, kp = kpf mod P. 

Finally, we will combine these results to deal with composite values (i.e. products 
of admissible primaries). 


5.1 N = P^, a power of an odd splitting prime 


We consider the case N = P^, with P an odd prime number, = 1, and /3 > 1. 

If the vector u = ^ ^ ^ generates a free submodule A in Z^, the condition 

(n n m n N) = 1 implies that n or m is coprime with P, and thus with N. In the 
case (n n N) = 1, n admits an inverse n~^ modulo N, such that n~^v is also a 
generator of A. We thus have the following; 

Lemma 2 For a primary value N = any free principal submodule A ofzj^ can 
be generated by a vector of the form 




(31) 


Then A = {nv, n G Z^r} is an eigenmodule of S (which we will note A = (v, A), 
with A the corresponding eigenvalue) iff the coefficient k (resp. k') is a solution of 

bk‘^ + (a — d)k — c = 0 mod N (32) 

(resp.) + {d — a)k' — 6 = 0 mod N. 


Generally, 6c ^ 0 mod P, so both 6 and c are invertible modulo N] we can then 
solve (^) provided D has a square-root in Zw, using the same formulae (^) as 
when working on K. Since P splits, there exist two such roots modulo N = P^, for 
any (3 > 1, which we note Pxps (cf. lemma |I|). Each of the above equations has 
then two solutions k± = 6“^{((i — a)/2 ± xat}, A:(i_ = c“^{(a — d)/2 ± xjv}, associated 
to the eigenvalues A± = (a -|- d)/2 ± xn- These solutions are of course redundant, 
since k±k'^ = 1 mod N. By convention, we can choose to generate 


the eigenmodule associated with A+ by the eigenvector 


A_ 


1 

k+ 

k'_ 

1 


(33) 


In the ‘caustic case’ be = 0 mod P, a square-root is easily extracted modulo P, 


D = 


a P d 


- 1 = 


a — d 


mod P 


(34) 


1 a — d , ^ 

we can take xp = - mod P 


12 







{P r\ D) = I implies (P n xp) = 1. Therefore, each of the equations (^) admits 
at least one solution modulo P, yielding two eigenspaces Ap^±, which can still be 
generated as in (|^ : the eigenvalue A+ = o (resp. A_ = d) is associated with 
kp = cj{a — d) (resp. k'_ = b/{d — a)). The extension from the case /3 = 1 to any 
/3 > 1 is done throngh the lemma ||. 

Although it is very general, the choice ( |^ ) for generators is not compulsory; 
for instance, if b is coprime to N, it is possible, and sometimes more convenient, to 

generate both eigenspaces by a vector of type 

In Appendix B, we similarly classify the eigenmodnles when N = P^, first for 
P an odd prime divisor of D {P\D 4^ = 0), then in the special case P = 2 

|[Tl[| . The analysis is slightly more involved than in the present section, but it only 
concerns a finite number of primes, which makes it much less ‘generic’ in the large-A^ 
limit. 

From the knowledge of free eigenmodnles for primary values N, one easily deals 
with composite valnes. 



5.2 N composite 

We treat the case N = 0^=1 = YYi=i ''^Ah Pi prime. If there exists a free 

eigenmodule A^r = {vn-,^n) modnlo N, then the submodnles Aj = {vi,Xi) modnlo 
Ni, given by 

Aj = Aw mod Wj Vi = l,---,n (35) 

it- — -iTat ” ” ” ” 

— vpf 


are free and invariant, hence every factor W is admissible. The Chinese Remainder 
Theorem |^, H] yields the inverse property: given free eigenmodnles Aj = (hj,Aj) 
modnlo Ni, the eqnations (^5[ ) have a uniqne solntion (Tat, Aw) modnlo N, cor¬ 
responding to a free eigenmodnle Aw, which is independent of the generators we 
choose for the Aj’s; indeed. Aw (although not Tw) is unchanged if we replace Vi by 
UiVi, where Ui is any invertible element of Zw^. 

( k' \ 

If both (6 n N) > 1 and (c n N) > 1, the generator ~ ) cannot be 


written in one of the forms (|3l|). However, the coordinates k'j^,k]^ can be chosen 
mutnally prime: indeed, the free nature of the submodule Aw is equivalent to 
(A:)y n A:w n A^) = 1. Therefore, G = (/cw H k'pf) is invertible in Zw, and Aw admits 
as generator vn/G, the coordinates of which are coprime. 

On the contrary, if (6 n N) = 1 (resp. (c n N) = 1), then, as in the previous 
section, k'j^ (resp. few) can be taken eqnal to 1 . 


Summary : Which N are admissible ? Let us now summarize our classi¬ 
fication (including the resnlts of Appendix B). If we decompose the discriminant of 
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S as in |24], D = Dgi, where the second factor is square-free, the following results 
hold: 

-for any prime P s.t. = 1) is admissible for any (3 (in the case P = 2, 

the splitting condition means D^i = 1 mod 8). We notice that = 1 is implied 

by = 1, but the converse is false. 


-for other primes, P^ is admissible if =0 and (3 is small enough; these cases 
only represent a finite set of primes (see Appendix B). 

-a composite number N = Pi'' is admissible iff each P^'" is. 

The number ctv of invariant strings for an admissible primary N = P^ depends 
upon (tnP). In the simplest case, = 1, cat = 2 for any (3 > 0. For a composite 
N = nr=i we have cn = nr=i CNi- 


How frequent are the admissible numbers among all integers, in the semi-classical 
limit ? According to Hardy and Wright |jl^ , a large number N is typically the prod¬ 
uct of few large primes (on the average, N contains oj{N) ~ log log iV prime factors). 
To be admissible, such a large N has to factorize into exclusively splitting primes, 
which represent statistically half of all primes. Therefore, admissible numbers be¬ 
come scarce in this limit. 


To illustrate our calculations, we represent some linear invariant strings for the 
transformation Sgg = ^ ^ ^ ^ 


3 2 


which is the simplest hyperbolic symplectic matrix 
in ri ,2 (cf @). Equations (|^) then reduce to = 3 or 3A:'^ = 1 modulo N. Since 


b = 1, all eigenvectors will be written in the form u = ^ ^ J . For the discriminant 

D = 3, the lowest splitting prime is 11, but 2 and 3 are admissible, although 
their powers are not. In figure 2, we show three invariant sublattices Aj^/N on T^, 
respectively for N = 11, 33, 121. These sublattices look well-distributed over as 
N increases. This remark leads us to the next section, where we investigate some 
semi-classical properties of invariant strings. 


6 Asymptotic properties of invariant strings 


We are interested in the way an S-invariant sublattice A^/N equidistributes over 
T^, for a large admissible N. To start with a simple case, let us assume that a 

generator u of A^v can be written u = ( ^ |, where A: is a solution of (|^ ) (this is 


possible, for instance, if (bnN) = 1). We consider a smooth observable f{x) on T^, 
and estimate its average over this sublattice. We are thus led to the sum 


N 


N 


< f >An= H f{nv/N) = ^ ^ ^ /p e 


227r- 


n=l 


N 


n=l 


(36) 
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by Fourier transforming the biperiodic observable. Permuting the sums, we see that 


1 

N 


N 

E 2 - 

e N 


n=l 


1 if p.v = 0 mod N 
0 otherwise. 


(37) 


Combining the above condition of ‘constructive interference’ with eq.(32), we obtain 


bpl — cp 2 + {d — a)pip 2 = Q{p) = 0 mod N. (38) 

For p ^ 0, this equation cannot hold in Z (equivalently, the isotropic set Q~^({0}) 
only intersects Z^ at the origin); indeed, multiplying Q{p) by 6 leads to 


bQ{p) = {Pib + - pID 


(39) 


and we know that D is not a perfect square in Q (due to the hyperbolicity of S). 
Therefore the equation ( |3^) implies p = 0 ov \Q{p)\ > N. If the largest eigenvalue of 
the quadratic form Q is noted qs = max | (6 — c/2)^ ± a/(6 — c/2)^ + D \, the second 
condition cannot be satisfied when p is inside the disk of radius y/N/qs- We then 
obtain the estimate 

|</>A. -/ol< E \fp\- (40) 

Equation (^) is meaningful if the Fourier transform of / is in Zi(Z^). However, 
since any continuous function / £ C(T^) can be uniformly approximated (i.e., for the 
II . Iloo norm) by trigonometric polynomials (see |l^, p.91), we obtain the following 
result; 


Proposition 3 For any infinite sequence of S-invariant sublattices \Am, Nadmissible], 
V/ G C(T2), hmAT^oo < / >A^= /o- 

The rate of convergence of this limit depends on the smoothness properties of 
/. Indeed, if / £ C“(T^) (i.e. f is a times differentiable), its Fourier coefficients are 
bounded as: 

\fp\<-^ for|^>l (41) 

where Ca depends on /. If a > 3, we obtain 

/ e C“(t 2) ^ |< / >^^ -/gl < . (42) 


When the generator v cannot be written in the form (31), the same type of esti¬ 
mates still hold. We outline the few changes needed for this case. By construction 
(see last section), N factorizes into two coprime numbers N = M M', such that the 
eigenmodule A^r admits a generator v verifying 


V = 



mod M, and 


V = 



mod M'. 


(43) 
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When estimating < / >An j the orthogonality condition p.v = 0 modulo N leads to: 


Q{p) = 0 mod M 
Q{p) = 0 mod M' 


Q{p) = 0 mod N. 


(44) 


The property Q{p) = 0 p = 0 still holds, so we obtain the same estimates as in 
eq.(@,||l). 


Taking the average of / over the lattices A^/N actually amounts to considering 
a set of atomic probability measures pN on T^; in this point of view, Proposition 
m means that such a sequence of measures converges, in the weak-* sense, to the 
Liouville measure on the torus, as N ^ oo, N admissible. We are interested in 
these atomic measures because they more or less correspond to the zero patterns 
of the quantum eigenstates we plan to build in the next two sections. Proposition 
m shows that the zeros of such eigenstates are uniformly distributed over in the 
semi-classical limit, which constitutes the main result of this article. 


Remark We notice that the above uniformity estimates were obtained quite sim¬ 
ply, and for all admissible N. On the contrary, the authors of |^| had to restrict N 
to splitting primes, and the uniformity property they obtained for atomic measures 
used non-trivial estimates of Kloosterman sums. The reason for this discrepancy is 
that these authors were considering families of atomic measures supported by indi¬ 
vidual periodic orbits modulo N, whereas ours are supported by whole eigenmodules, 
which are unions of ‘ideal’ periodic orbits pi]. 


7 Prom classical invariant strings to quantum 
eigenstates 

Now that we have studied some asymptotic properties of invariant strings, we choose 
one of these strings for an admissible value N and express the quantum eigenstate 
with corresponding (lattice-like) zero pattern as a pure theta state in some invariant 
space 'Hn,k- We then do the same for another type of invariant sublattices, namely 
square sublattices, which appear when iV is a perfect square. In the subsequent 
section, we will then build a larger family of eigenstates for the same N, allowing 
their zero patterns to be crystalline instead of lattice-like. 

7.1 Eigenzeros forming a lattice, N admissible 

We start from an eigenmodule Aiin = (u, A), obtained in section (P) . We can then 
search a corresponding affine invariant lattice A, i.e. find a real vector Vq s.t. A = 
Vo + Alin is itself invariant through S. We will not make here an exhaustive list 
of such affine invariant strings, but rather restrict the entries of Vq to be integer 
multiples of V/2. Each eigenstate we will exhibit will then belong to one of the four 
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spaces TCn,k, with k =0, 0, 1/2, i/2. We divide the study according to the parity of 
N, for reasons to become obvious. 

Let us start by considering N odd and admissible; we select a free eigenmodule 

' k' ' 


Alin with a generating vector v = 


k 


where k and k' are coprime. Since N 


is odd, we may choose the entries /c, k' to be of opposite parity. Indeed, suppose 
they are both odd; the following transformations yield a new generator of Ann in 
the desired form, 


V = 


k 


k' + N 
k 


^ 1 


k' + N 
k 


where G = {kr\{k' + N)). 

(45) 


In that case, we can find two integers u, v, s.t. the matrix 

S,-t = 


u k' 

V k 


(46) 


belongs to ri ^2 (see eq. 
( 0 ' 


- 


1 


. This matrix maps the vertical string Aqo generated by 
modulo N, onto Ann. Since Aqo is associated to the antiperiodic basis 


state |^{Ar-i )/2 = l/2)Ar6’ eigenstate associated with Aim is itself antiperiodic, 
and reads as 

\^)n,6 = Us,\1/2)^^^ . (47) 

The zero pattern of this state consists in the lattice Aim/A'v^, obtained from Aim 
through a dilation and reflection w.r.to the real axis (or complex conjugation). Since 
N is odd, the four strings 


A 


lin? 


A^ — Alin + 


N/2 

N/2 


Ai — Aim + 


N/2 

0 


— i^lin + 


0 

N/2 


(48) 


are all different. Since S G ri _25 A' is invariant, and yields a periodic eigenstate 

\'‘P')n,0 = UsJQo = 0)Ar,0 • (49) 

If moreover S is congruent to the identity matrix modulo 2 (equivalently, 5 G r 2 , 
the level-2 principal congruence subgroup [^]) , then Ai and Ai are invariant as 
well, and yield respectively the eigenstates 


\'>Pi)n, 1/2 — Us^\qo — 0)Ar,l/2; 

l^i)Af,i/2 = t^Ss|l/2)Ar,i/2- 


(50) 


When N is even and admissible, the coordinates k', k of the generator v of Aim 
can still be chosen coprime, but we cannot modify their parity through the trick 
of eq. (|45|) , so the matrix of equation (|4^) is not always in ri_ 2 - Besides, the 


four sublattices given in eq.(^) are not all different. For example, if k, k' are odd. 
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then Alin = A', Ai = Ai are invariant and the corresponding eigenstates are both 
antiperiodic, 


N,d ~ ^Ss\'^/‘^)n,1/2^ (51) 

= C^5s|0)Ar,1/2- 

Similar phenomena occur for different parities of k, k'. For instance, if k' is even, 
then Alin = Ai, A' = Ai are both invariant, and the two eigenstates are in 'Hn^ii 2 - 

To summarize, in the case where N is odd admissible, we can build from any 
invariant linear string A^n both a periodic and an anti-periodic eigenstate (plus two 
other eigenstates iff S' G r 2 ); when N is even, we can build from Aun two eigenstates 
which belong to the same space with the Floquet parameter k depending on 

the invariant string Aun we consider. In hgure 3, we draw the Husimi functions of 
two lattice-like periodic eigenstates of for the two admissible values N = 33 

and N = 121; these eigenstates are built using the invariant strings drawn on figure 
2 (center and right). 

The following section deals with a slightly different type of eigenstates in the 
particular case when A^ is a perfect square. 


7.2 Eigenzeros on a square lattice 

Irrespectively of the number theory leading to free eigenmodules, a simple integer 
lattice shows up when A^ is a perfect square, say N = M^: namely, the square 
lattice Agq = M{Z iZ), invariant through any modular transformation S, and 
containing N points modulo N (Agq is also a submodule of Z^, but of type 2 
and not free). Since the unique state of TC^ g has a simple zero at the origin, the 
function z i—> (Mz|l/2)^g has zeros on the square lattice of side l/M\/2, i-e. on 
Asq/N\/2. A straightforward calculation shows that this function is the Bargmann 
representation of a state {z\ijjsq)N,K, periodic if N is even, anti-periodic if N is odd. 
Writing down the theta series, we easily get the decomposition of \’tpsq)N,K in the 
position basis, 

1 1 /2 

= 75 S <“> 

with K = 0 (resp. k = 0) for N even (resp. odd). These states are not pure 
theta states according to the dehnition we gave in section ^ since they are not 
built through the action of Mp(2,Z) onto some basis state of 'Hn,k- However, they 
transform simply, as shown below. Taking any transformation S G ri^ 2 ) we label 
the operator Us according to the space it acts upon: 

(z|[/J'^-|^gq)jv,. = {Mz\Us^’°\l/2)^^q = e2--6(Mz|l/2),^g (53) 


with K = 0 or K = 0 according to the parity of N. The first equality is due to 
the linearity of the transformations in eq.(^), the second is obvious is a 
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one-dimensional space on which Us is an endomorphism). The eigenangle Ug only 
depends on the transformation S. 

For the same N = we can build a second eigenstate associated to the 

' M/2 ' 


invariant lattice A'q = Agq -|- 
it is periodic for all M : 


M/2 


The quantum state can be built as above, 


M-l 


IV^sq)Af,0 /VT ^ ' |^/Af)jV,0 • 




(54) 


Proposition 4 For any S G ri ^27 for any perfect square N, the states I'i/sq) n,k, iV'sq)A^.o 
are eigenstates ofUs- 

In figure 4, we draw the Husimi function of the periodic eigenstate |'!/sq)i 2 i,o for 
N = 121. 


8 Crystalline eigenfunctions 

In the last section, we built eigenstates for which the zero patterns were classically- 
invariant lattices A, which filled the unit torus uniformly in the semi-classical limit. 
We would like to find a larger class of eigenstates with the same property. 

The authors of derived semi-classical equidistribution properties concerning 
the Wigner distribution of eomplete orthonormal families of eigenstates, although 
they had to restrict N to splitting prime values to obtain these results. Equivalently, 
they proved that for these values of N, Schnirelman’s theorem is valid for any 
sequence of eigenstates; their proof relies on the bilinearity of the semi-classical 
measures (Wigner of Husimi) w.r.to the quantum states. 

On the opposite, we are not directly concerned with the properties of the Husimi 
densities, but rather with the atomic measures associated to their zero patterns. 
Since these patterns are very non-linear objects, we are not able to generalize the 
above uniformity properties to complete families of eigenstates. However, we can 
obtain a partial generalization from last section: when N is not a prime, it is 
sometimes possible to build eigenstates with zero patterns forming erystals instead 
of lattices, whose periodicities correspond to non-free eigenmodules of Z^. 

These crystalline states can be built through different ways, depending on the 
value of N. 

8.1 N admissible 

To simplify the presentation, we will give details only in the case where N is odd 
admissible, and moreover we first consider an invariant linear string A = Aun = 
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with eigenvalue A = a + bk. Therefore, the 


{v, A) such that v is of the form 


k 


matrix of eq. (|4q) can be written 


S^ = 


0 1 

-1 k 


(55) 


We note the vertical integer translations of A (which are, in general, not invariant 
under S): 


Ai = A + 


0 

I 




Aoo + 


-I 

0 


for I G Ztv- 


(56) 


By linearity, S maps the string A; into a parallel string A;/. More precisely. A;/ = 
Ayi, where X' = d — bk is the inverse of A modulo N. To these translated strings 
are associated the anti-periodic pure theta states: 


|V>z),v,o = ^s,|l/2-//iV)^,o ■ (57) 

These states are not eigenvectors of Us, but they form a new orthonormal basis of 
on which Us acts quite simply: 




N,0 


(58) 


The phases (t{N, I) can be derived by writing down explicitly the matrix elements 
of Us in this new basis. This leads to more or less complicated formulas, some of 
which are given in Appendix A (see also j^, 20, for the case k = 0). 


To build a basis of eigenstates, we use the partition of Zy into disjoint cosets 
(or orbits) Oj = {A'^ lj,s = 0, - ■■ ,pj — 1}: Ij is an arbitrary point on the orbit Oj, 
and pj the period of the orbit. Using the usual notation ())(n) for Euler’s totient 
function |jT^, we have necessarily pj\(j){Ni.) (we recall that Ni = N/{N Hi)). 


Starting from an element L of the coset Oj , we can build a family of pj orthonor¬ 


mal eigenstates of Us- According to eq. (|58[) , IV’ij )y 5 and its successive images under 
Us are all eigenstates of U^\ with the same eigenvalue exp( 2 i 7 r(Tj), where 


cr,- 


Pj-l 


s=0 


(59) 


We therefore obtain pj anti-periodic eigenstates |</>j,r) 7 vo with eigenvalues 

^2m(r+aj)/pj ^ for r = 0 , • • • , Pj _ 1 : 




1 

vfo s =0 




(60) 

with 

Vr = 0, 

arbitrary, and 


< 

II 

• • • u ■ - 2 

■:P3 fo-l -1 

_ ^-2iTT{r+aj)/pj ^2m(i{Nlo) ^(r) ^ 



20 






To each coset Oj thus corresponds a set of pj orthonormal eigenstates, with the 
eigenvalues written above. Two eigenstates built from different cosets are orthog¬ 
onal, since they involve disjoint sets of orthonormal states |V'z)Aro- However, they 
can share the same eigenvalue ||^, |^. By completing the above procedure for all 
cosets Oj (including the trivial one Oq = {0}), we obtain an orthonormal basis of 
(possibly degenerate) eigenstates. 

In general, these eigenstates are not pure theta states, and their zero patterns are 
not known analytically. However, if N is not a prime and (IjCiN) > 1 (this gcd does 
not depend on the element Ij we have selected in Oj), then the Bargmann functions 
of the eigenstates \4>j,r), ^’ = 0, • • • — 1 inherit partial quasi-periodicity properties 

from their components (^IV’z)- To prove this, we generalize the fine periodicity of 
basis states (eq. (p^) to all pure theta states: 


( 2 + 


= 2 i 7 rgj 


(61) 


If we apply this formula to the components of \(f)j^r) and iterate it, we obtain 


/ , + i^l ; 


ii)N.d ~ h)N.d ( 


g-7r/i(l-|-ifc)\/22 


e 2 JV 


2i7r^A^'^ Ij 

e N 


(62) 

The only s-dependent factor is the last one. If we choose p = Ni., this dependence 
is removed, so that the eigenstates \(pj,r)]sfo themselves share the above periodicity 
property, 


/ , Ar + \ / IJ, \ / -TTNi.{l+ik)V2z 7rAr2 (l+fc2)/2Ar 

+ \^h,r)N,0 = {z\^h,r)N,0 ^ e h 

(63) 

As a consequence, the N zeros of {z\4>j^r )q form a crystal, whose lattice is generated 


by the vector (in complex representation) (— l-|-iA:)/(A^nZj)\/2, and contains (NCilj) 
points; each fundamental cell contains Ni^ zeros, but we have no information about 
their locations within the cell. 

If we rename our basis of eigenstates as 




(64) 


we obtain a basis of orthonormal eigenstates / = 0, • • •, — 1} s.t. each 

zero pattern of is a crystal generated by (—1 -|- ik)/{N n l)V2. 

Note that if we perform the same construction of eigenstates using a different 
invariant string, we will a priori obtain a different orthonormal eigenbasis; this is 
not inconsistent, considering the frequent degeneracies of the spectrum [^]. There 
can actually be redundancies between two eigenstate bases built from two different 
eigenmodules A^^). Indeed, if, for a certain divisor g of N, we have = 
mod N/g, then the crystals of the eigenstates with g\li, g\l 2 , are 

both supported by the same lattice; if these eigenstates share the same eigenvalue, 
nothing prevents them from being rigorously identical. 
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To deal with an invariant sublattice with generator v = 



i.e. the case 


where both (b n N) > 1 and (c n N) > 1), we cannot restrict ourselves to vertical 
translations of A, as was done so far. To build a complete eigenstate basis, we have 
to consider all possible translations of A of the form Vq + A, with Vq an integer 
vector. Thus, Vq should take values in a fundamental cell of the sublattice A, i.e. a 
representative of the quotient (NZ)‘^/A. Such a fundamental cell can be obtained by 
the same procedure as above, i.e. by 5ij-transforming (using eq. (^)) the vertical 

strings Aqo + ( ^ ),for/ = 0,---,A' — 1. The construction of eigenstates then 


0 


proceeds as above, using the inverse eigenvalue A' = k{au + bv) — k'{cu + dv). 

This construction can be generalized straightforwardly to obtain a basis of pe¬ 
riodic {k = 0) eigenstates for the same odd N, by using integer translations of the 
string A' (see section O). One obtains a basis of orthonormal eigenstates |<h;)Ar,o 
similar to the one obtained above (eq. (|^ ). 

As an illustration, figure 5 shows the Husimi functions of two crystal eigenstates 
of lAsfig for = 33 = 3 X 11 (both 3 and 11 are admissible). On the left, the zero 
pattern is periodic w.r. to a lattice of order 11 , whereas on the right, the periodicity 
lattice is of order 3. 

To deal with an even admissible N, one just needs to take care of the Floquet 


parameter k the same way as in section 7.1 


Let us turn back to the orthonormal basis of antiperiodic eigenstates {1$^)^ q, I G 
I^n}- Noting g = (NCil), the state 1$;)^ 5 has a zero pattern of period {—1+ik)/g\/2, 
which generates an invariant string of g points on Tc. Through a dilation of factor 
gV^, this linear string is mapped into a free principal eigenmodule of Averaging 
an observable / over the crystal of zeros, we can therefore apply equidistribution 
estimates identical to eq. 


l|,12), after replacing by s' in the formulas. The zero 
patterns of are thus well distributed over Tc when g is large. Unfortu¬ 

nately, the eigenstates featuring such regular patterns are not very numerous: for 
any divisor g of N, the number of integers I G Zjv s.t. {N D 1) = g is given by 
(p{N/g), and Euler’s totient function varies like (/>(n) ~ ^nj'K^ on average [12|. As 


a matter of fact, the eigenstates |‘hz)jv 6 no a priori crystal periodicity (i.e, N 
and I coprime) represent a non-negligible fraction <j){N)/N of the whole basis. 


Finally, we notice that the quasi-periodicity of the Bargmann eigenfunction (^) 
is mapped to an exact periodicity of the corresponding Husimi density. Therefore, 
the uniformity estimates (^,^) apply as well if one averages / over the Husimi 


density itself. From the results above, we therefore deduce the following proposition: 


Proposition 5 Consider an infinite sequence M = {N G N, N odd admissible} 
and a corresponding sequence of crystal eigenstates {iTjAr = | ^i)n, glATeA/" such that 
g = {N r\ 1) —>■ 00 . Then both the Husimi measures and the atomic 

measures {gN}NeJV describing the zero patterns of H^j^ converge weak-* to the 
Liouville measure on as N ^ 00 , N G M. 
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8.2 N a perfect square 

We consider again the case N = and construct an orthonormal family of 

periodic (k = 0) eigenstates. In section |7.2|, we built the eigenstate associated with 
the lattice A'q = M(Z + iZ) + M/2(l + i), 

(^IV’sq)Af,o = {Mz\qo = 0)i,o . (65) 


We generalize the above formula to include states with zeros on translates of Ag^/N^/2. 
One easily checks that the state \'iP)n defined by {z\'iP)n = {Mz\qo = K 2 )i,k belongs 
to Ti.N,o iff ^ is of the form k = (ni + m 2 )/M with ni,n 2 £ its zero pattern 


is associated to A' + 


n2 

ni 


This family of periodic states is globally invariant 


under Us, according to a generalization of eq.(53), 


{Mz\u'^^'‘ 


1^2/1,K — e 


2i7rc7« 


{Mz\k2] 


1,k' 


with K = S K mod [1, i]. 


( 66 ) 


Moreover, these states are mutually orthogonal (this can be seen by decomposing 
them on both position and momentum bases in Ti-Nfi), so they form an orthonormal 
basis of U-Nfi- Starting from a certain of the above type, member of the coset 
Oj = {S'* i = 0, - ■■ ,pj — 1}, we can then build pj periodic orthonormal 

eigenstates \4>j,r)N,o of Us, involving linear combinations of (Mz|k 2 ■ ^s in 
last section, we finally obtain an orthonormal basis of eigenstates. 

If the coordinates of have the gcd (ni n 77-2 n M) = g, then all the states 
I^ 2 ^)ik (0 share the following quasi-periodicity relations; 


{M{z + l/gV2)\K^^)^^^(i) 
{M{z+ i/gV2)\K^2^)^^^(i) 


^^{Mlg)^+nV2z{MUg) (^M z\4'>) (67) 

^^{Mlg)^-izrV2z{MUg) . 


Since these relations are (7)-independent, they apply to each eigenstate |0j,r)w,o as 
well, the N zeros of which belong to a crystal supported by the square sublattice 
in Tc of side ljg^/2 (each fundamental cell of the crystal contains {M/g)‘^ = N/g‘^ 
zeros). 

In figure 6 (left) we plot the Husimi function of an eigenstate of Usf^^ for the 
square value N = 81. The starting Floquet parameter is = (3 — 3i)/9, so we 
obtain a square crystal of side length 1/3. 

In the case N = odd, a similar construction leads to anti-periodic eigenstates 
with the same type of zero patterns. 

As in the former subsection, the crystal of zeros will be well-distributed on the 
torus when g is large, which happens for relatively few eigenstates of the basis 
(namely, when both ni and re 2 are multiples of g). 


8.3 N with an admissible divisor 

Since we are in search of families of eigenstates displaying semi-classical (i.e. N 
oo) properties, the admissibility of N seems a very restrictive condition. Indeed, 
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this property implies that all prime divisors of N are admissible themselves, which 
represents asymptotically half of all primes. Therefore, the admissible numbers form 
a set of vanishing density among all integers in the semi-classical limit. For this 
reason, we now use the ideas of the former sections to build crystalline eigenstates 
in the case where N is not admissible itself, but has an admissible divisor M, say 
N = M m {M is supposed maximal). We will build eigenstates at this value N^ 
for which the zero pattern verifies the same periodicities as eigenstates built at the 
value M (we restrict ourselves to odd N and anti-periodic {k = 0) states, which are 
the most amenable cases). 


The construction still uses classical results. We start from a free eigenmodule 


A modulo M, generated by v 


1 

k 


(we choose this form for v to simplify 


notations, so that the matrix can be written as in equation (pSD). Obviously, 
the lattice mA is invariant modulo N (it is a non-free principal eigenmodule of S in 
Z^). We now consider a free module A° in Z^, which contains mA as a submodule 
(for instance, we can take the string generated by the same v modulo N). A° is 
not invariant under S, but its successive images A"" = 5” A° all contain mA. We 
note prot the least integer s.t. = A° (prot is typically of order m). This 

period corresponds to a cycle of ‘rotations’ of A° under S, as opposed to the cycles 
of translations considered in the last sections. 


Quantum-mechanically, the state 




( 68 ) 


and its successive images Ug\'tp°)^Q are pure theta states associated respectively 
to the lattices A”. As a consequence, each of the states Us'^\'ip °)n q has M of its 
zeros on the string A/M^/2. Besides, |'0°)Ar 5 ~ phase a^ot 

depending on S and A°. We can then build a family of prot orthogonal eigenstates 
with eigenangles for r = 0, • • •, (prot — 1): 


Plot -1 


\ 4 >- 


r/N,0 


E 


n=0 


us-\r)M.6 • 


(69) 


This equation is similar to the formulas (58,60) we derived in the case of N ad¬ 


missible. However, as opposed to section 8.1, these orthogonal eigenstates are not 


normalized., since the components t^5lV’°)7vo mutually orthogonal; some of 

the |</>r)jvo even be null vectors, if |V’‘’)Aro orthogonal to the corresponding 
eigenspaces (we expect this situation to be non-generic). 

The Bargmann functions of these pmt eigenstates all vanish on A/M-v/2, since 
all their components do. More precisely, these components (and thus all the eigen¬ 
states) possess the same quasi-periodicity w.r.to this sublattice: 


Vn, {z + m 


— 1 + ik 
Nt/2 


\Us^\r)N,6 = i^\Usir)N.6 i-irem 


m p—-Km{l+ik)\/2z 


(70) 
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We have therefore built a family of orthogonal eigenstates whose zero patterns 
are crystals generated by For the inverse Planck’s constant M, there was a 

unique state with such a property. In figure 6 (right), we plot the Husimi function 
of a periodic eigenstate of Usf^^ in the case = 55 (11 is admissible, 5 is not); the 
invariant string A modulo 11 is shown in figure 2 (left). 


Further on, we can build eigenstates with coarser quasi-periodicities, by compos¬ 
ing ‘rotations’ and translations of A°. Classically, the integer string A° ^ ^ = 

A° -|- loVoo modulo N is mapped under successive actions of S into sublattices of the 
form AP -\- Vj, where Vj are integer vectors. We know that will be a translate of 
A° iff j is a multiple of prot ■ 

We can then use the results of section (^) to describe the cycle of translates of 


A° under the action of 5rot = 5'^“'°*; this matrix maps A° -|- IqVoo into A° + loVoo, 
where is an eigenvalue of S'rot modulo N. Without deriving explicitly the 
matrix S'rot; we know that A(.o^- = mod M, where X' = d — bk describes the 
translates of A modulo M (cf. eq. (|^). We can therefore decompose into 
cosets Oj = {A'rot/j, s = 0,... — 1} (including the trivial coset Oq = {0}). For 

each coset, the total period of the cycle {S'” {A° + /jUoo)} is thus ProtPj- 

Considering the coset Oj, the pure theta state we start with is = 


UsJl/2 — Ij)Q- Its images through Us are associated to the lattices A^ -|- Vj, 
so they have a priori no common zero with However, as in section (8.1), 

these states verify some non-trivial quasi-periodicity properties, cf. eq. 10), 




m p-7rm(l+ifc)C22: p-2i7r- 


(71) 

and the only n-dependence appears in the last factor. From this sequence, PmtPj 
orthogonal eigenstates can be built using formula (^9[), with initial state ) 7 v 5 
period prot Pj (tbe phase a now also depends on the coset Oj). A priori, some of these 
states could be null, for the same non-orthogonality reasons as above. However, if 
such an eigenstate is not null and g = (Ij n M) > 1 , its zero pattern is a crystal 
generated by (—1 -|- \k)lg\f2. 

For each coset Oj , we have obtained a family of pj prot orthogonal eigenstates 
(some of which can be null), with zeros on a crystal generated by (- 1 - 1 - ik)/{lj H 
M)\/2. These crystals are based on the same lattices which appeared for the value 
M, but now each fundamental cell contains m{lj n M) zeros. Unfortunately, the 
structure of the whole family of eigenstates is less clear than in section | 8 . 1 | : on the 
one hand, we do not know the norms of the eigenstates, even their very existence; 
on the other hand, two degenerate eigenstates built from different cosets are a 
priori not mutually orthogonal. Actually, this lack of information is not surprising, 
since we have built prot times too many eigenstates through our procedure, so that 
dependence relations between them must exist. However, as in the last section, 
equation dZH) ensures that Proposition ^ applies to sequences of non-null eigenstates 
5 for which g = {Ij n M) ^ oo. 
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9 Concluding remarks 

Both the linearity of the transformation S and the classical invariance of certain 
sublattices (strings) of transpose to quantum mechanics, when expressed in the 
Bargmann-Husimi representation. On the time-dependent side, the quantum map 
Us acts classically on string-like zero patterns (i.e., zeros of pure theta states); on the 
time-independent side, it admits families of eigenstates for which the zero patterns 
are crystals supported by classically invariant sublattices. 

The lattices supporting these crystals (thus the crystals themselves) equidis- 
tribute on the torus when the number of lattice sites increases to infinity, and the 
equidistribution of the corresponding Husimi densities then results as a byproduct. 

However, crystalline eigenstates do not emerge for all values of N, and when 
they do, they concern only isolated states, whereas eigenspaces are very often multi¬ 
dimensional, due to arithmetical degeneracies of Us- Obviously, a linear combina¬ 
tion of two eigenstates with crystals supported by different sublattices will not be 
crystalline, nor will their zeros be a priori well-distributed on Tp: the equidistri¬ 
bution of the zeros is not a linear property. Therefore, the eigenstates for which 
we have a nice zero pattern (i.e., a well-distributed crystal) only represent a few 
elements of the whole set of eigenstates; they do not even form a basis (indeed, they 
usually span a rather small subspace). Moreover, such crystal features certainly 
disappear as soon as one considers a non-linear perturbation of S. 

In reverse, for many other eigenstates of Us that we have computed numerically, 
zeros show no particular periodicity, yet remain more or less well distributed on 
the torus. Surprisingly, this property persists even when the Husimi function of 
the eigenstate has a ‘scar’ on a classical fixed point (see figure 7). Therefore, the 
equidistribution property seems to be even more general for the zeros than for the 
Husimi density itself, but a proof will probably involve a different approach from 
that to Schnirelman’s theorem. 
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also indebted to S. de Bievre, M. Degli Esposti, S. Graffi, P. Leboeuf and S. Zelditch 
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Appendices 


A Derivation of the matrices Us 


In this appendix, we derive the entries of the N x N matrix which represents the 
operator f/5 mapping the basis {|gj)Ar,K, j = 1} to the basis / = 

0, • • •, V — 1} (see eq. (|^). The Floquet parameters k, k' G [0,1)^ are related by 




(72) 


where the integers no,mo are called winding numbers We will only consider 
hyperbolic matrices S, so the entry b cannot vanish. The formulas obtained are 
generally used when Us is an endomorphism, i.e. n = k! , but we do not need this 
property in our calculations. For the periodic case (k = k' = 0), the matrix elements 
were derived by Berry and Hannay |^] for any N. We recall how to obtain these 
elements for general Floquet parameters, and write explicit formulas in two cases, 
first when N is coprime to b (cf. our example S'fig), then when N divides b (which 
is useful for section O). 

The Berry-Hannay approach consists in extending the action of Us from square- 
integrable functions to tempered quasi-periodic distributions 5'(IR) (see also |^|). 
Precisely, from the image of a unique Dirac peak. 


Us\q) = I dq' {q'\Us\q) W) (73) 

given by the kernel ([T^), one derives the image of the quasi-periodic distribution 
\qj)N,K- From equation (^), we decompose the translation index as m = bu + k, with 
n G Z, k = 0, ■ ■ ■ ,b — 1, and obtain, using Poisson’s summation formula. 


Us\qj)N,K = 


1 

ViNb 


b-l 

nez A:=0 


aqj+bni/N+n/N). 


(74) 

In this equation, we already see that the image of \qj)N,K is a sum of regularly 
spaced Dirac peaks. We now show that this image is itself quasi-periodic. For this 
matter, we compare the coefficients in front of {q^ + n/N) and \q'^ + n/N + 1), noting 
q'o = + bni/N. 

To avoid too cumbersome notations, we suppose from now on that {b H N) = 1. 
The sum over k then yields Q 


Us\qj 


•N,k = 


Cn 

ViN 




YfaNiaN\bfn^ ^i^{dq'^-2qjq'+aq^) |^, ^ ^ 


where Cjq is a number-theoretical phase (see eq.^), and {aN\b) is the integer 
inverse of aN modulo b (this inverse exists and is unique since {aN Cib) = 1). The 
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exponents in front of + n/N + 1) and + n/N) differ by 


-jaN{aN\bf{2nN + N^) + {2d{q'^ + n/N) + d - 2qj) 


(76) 


ITT 


= 2i7rK'^ + 2i7r(cj + mo) + ~ a,N‘^{aN\b)‘^). 

Finally, a simple argument shows that (d — aN'^{aN\b)‘^) /b is an even integer. 
Therefore, we obtain a finite sum of quasi-periodic distributions. 


Us\q,)N,. = -fe E* \q' = q'o + n/N)^^,,. (77) 

V n=0 


7V-1 


To be complete, we give the value of the global factor Cn, 
if 6 is odd, Cn = 


(78) 


if b is even, Cn = ( ^ ) 


where is the Jacobi symbol, i.e. an extension of the Legendre symbol of eq.(^) 
to non-prime b (see im for details). The unitarity of the matrix Us between the 
Hilbert spaces 7iN,K and Ti.N,K' can be proven by decomposing L 2 (R) into a direct 
integral of spaces Hn,k [H]: it then follows directly from the unitarity of Us acting 
on L 2 (K). 

Identifying each space 'Hn,k with , we now express the matrix elements be¬ 
tween the states N,K'{Qj'\ and \qj)N,K for j,/ = 0, • • • , A^— 1. Due to the assumption 
(5 n N) = 1, b admits a unique inverse b~^ = {b\N) modulo N. With this assump¬ 
tion, we can also choose the integer b~^ s.t. a b~^ and db~^ are even numbers. After 
a few calculations, eq. dH) yields the following matrix elements. 


N,^'{qj'\Us\qj)N,K = CArexpj^ j'‘^{db ^l2)-jjb ^+j‘^{ab ^/2) 


(79) 


+j'{ k[ + mo - db ^no)+j{b - ACi)j | 


with the prefactor 


Cn 


|-^((iK 2 ^ — 2k 2«^2 + ^^i) “ '—aN{aN\b)‘^n\ 


In the periodic case k = k' = 0, we recover an expression very similar in form to 
the one in the continuum, except for the prefactor (cf. ill)- Since we are mainly 
interested in Us being an endomorphism, we apply the above formula in the case 
K = k', which is possible for a finite number of parameters k, indexed by the winding 
numbers no,mo 


K2 

Kl 


2 — a — d 


d-l -b 

—c a — 1 


Ho 

mo 


(80) 








In the main text, we have only built eigenstates with periodicity angles k = 0,1/2, i/2, 0 = 
(1 + i)/2, but once given the matrix S, we can easily extend our results to all in¬ 
variant parameters Kn^^mo- 


From eq.(|75D on, the above calculations depended on the assumption (bON) = 1, 
which for instance is true for all values N in the case 6 = 1. If now (6 n N) > 1, the 
image of \qj)N,K is still a superposition of quasi-periodic distributions in but 

now some elements n ,k' {q'jWs\qj)N ,k vanish. We are led to consider the case A^|6 in 
section (0), in the course of deriving the matrix elements i\Us\'4^i) q- Indeed, 
group properties of Ug lead to 




(81) 


(the ± sign does not depend on the bra-kets considered). For simplicity, we take N 

odd, and suppose that can be written & = ( 

\ —1 K 

to eq. (|3^) , the matrix S = S writes 


with k even. According 


S = 


X’ = d-bk (5N 
—b X = a + bk 


(82) 


where A is the classical eigenvalue associated with ^ ~ ^ ^ J modulo N, X' its 

inverse, and P = [bk‘^ + {a — d)k — c]/N is an integer coprime to Ab By the same 
type of calculations as above, we obtain 


n,k' {Qj'\Ug\qj) — 


= ^exp |^[A'm2 + 2m{q,X! + Pki - qj,)] 

(83) 

Here, (.)^ means the average over the integer m. This average is non-vanishing iff 
qjX' -|- Pki — qji is integer, which corresponds, for each j, to a unique index j' = jg 
modulo N. Thus, the matrix acts, up to phase factors, as a permutation between 
both bases (cf eq.(|^)). Using the winding numbers 


S K = K + 


n 

m 


(84) 


the state \qj)N,K is mapped by into \qjs)N,K’: where jg = X'j + n. The phase 
is given by 

N,.'{qjs\Ug\q,)N,. = Ck 

— gi7r/3(AKf-2feftiK2) g-i7r6A'K|/Af P+2j{h+K2)] 


where the constant Ck = e if 6 (and thus /3) is odd, Ck = (y) ei7r(A'-i-i)/4 
if 6 is even. The formula simplifies greatly in the periodic case {k = k' = 0). 
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Since most of the eigenstates built in the text are antiperiodic {n = k' = 0), we 
now link the notations above to the ones in section e, i-e. we map the indices 


j,js to the indices I, I' of eq.(^), and adjust the sign of the image state: 
on the classical side on the quantum side 

^ ^ ■ ki)iv,6 = \j/^ + V2iV)iv,o = 11/2 - VA^)iv,o 

A'+/3-l , 


J = 


-I 


js = A'j + h, l' = X'l kis)Ar,6 = (-1) " |l/2-0Ar,( 


( 86 ) 


B Free eigenmodules, continued 

As a complement to section (0)> we classify the invariant strings modulo N, first 
when is a power of a prime divisor of D, then when is a power of 2. We use 
the results of section (5.1), in particular lemma ^ still applies. 


B.l N = with P an odd prime, P\D 

Let P > 3 be a prime divisor of P, a > 0 be the largest integer such that P°‘\D. 
We seek to diagonalize S modulo N = , for /3 > 1. 

We first assume that P does not divide b and c simultaneously, therefore all 
eigenmodules can be obtained by solving one of the eqs. (|3^ , so the problem reduces 
to extracting a square-root of D modulo N = as in section (5.1). 

For (3 < a, the roots are the elements x = X pK^+i)/^]^ where the brackets [ ] 
mean the integral part, X taking any integer value. 

For (3 > a, D = P°‘6 is a square in Zp/3 iff a is even and = 1. In those 

cases, its square-roots take the form x = Py P'^/‘^ + X P ^“'2 , where y verifies = 
6 mod P^“", and X can be any integer. 

In both cases, to each square-root x modulo P^ there corresponds one free 
eigenmodule A. 

In the case where P divides b and c simultaneously, we decompose S as 

5 = ± I + P^ I ® -f 1 = ±/2 + P^M, (87) 


( 88 ) 




with the condition M ^ 0 mod P. We obtain the following constraints on M: 

det5 = 1 (e + /i) = =F P'^{eh — gf) 

r /. , u\‘i 

(e + h) is even and D = P^"^ 




which implies 2'y < a. As long as [3 < 'j, every vector of is invariant modulo 
P^, so A can be any free submodule of type one. For larger values of /3, we need to 
diagonalize M modulo P^“'^, that is, to solve the remainders of eq. (32), 

(89) 


d-7 


or 


fP + (e — h)k — y = 0 mod P' 
gk'^ + (h — e)k' — / = 0 mod P^~'^. 
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This amounts to extracting a square-root of the discriminant of M, 5' = 

(this discriminant is related to the 5 defined above by 6' = 

If / or 5 is coprime with P, the solutions of (|^) will depend upon : if 


= 1, we apply the results of section (5.1) and diagonalize M for any /3; if P\6' 
(that is, a > 27 ), the first case considered in this appendix provides eigenvectors 
-at least up to /3 = a — 7 , 

-for all /3 iff a is even and = 1 ; 


If = —Ij there is no eigenvector for /3 > 7. 

The case f = g = 0 mod P is actually solved easily: since 7 is maximal, 
e = —h ^ 0 mod P, so the equations ( | 8 ^) have two solutions for any /? > 7, using 
lemma |I]. 

To conclude, we recall that any solution of (^) modulo pl^~'^ provides P'^ dis¬ 
tinct invariant strings modulo P^. 


B.2 N a power of 2 


The prime number 2 has to be distinguished from the other primes, especially 
when one considers quadratic equations | 11 ], which play an important role in the 
diagonalization of S. We separate the two possible forms of S modulo 2 . 


In the case S = 


S admits the unique eigenvector 


mod 2 . 


Since b, c are odd and a — d even, eigenvectors are obtained by the extraction of a 
square-root of D modulo 2^, for any /3 > 1. 

If (a-|-d)/2 is even, then D = 3 mod 4 is not a square |Tl|, so there is no solution. 

If (a -|- d)/2 is odd, we write D = 2“5, with 6 odd, a > 3. As long as (5 < a, we 
have the trivial square-roots x = X for any integer X. When /3 = a -|- 1, 

there are solutions iff a is even, x = 2“/^ -|- X for [3 = a + 2, the existence 

of solutions requires 5 = 1 mod 4, and for j3 > a + 3, the condition to extract a 
square-root is d = 1 mod 8. In the last two cases, all square roots are of the form 
X = ±2“/^y -|- X where ^ 6 mod 2^“". 


The alternative case 


S = ± 



+ 2 ''^ 


e / 
9 h 


= ±l2 + 2T'M 


(90) 


leads to similar phenomena as for odd primes. 

For /3 < 7 , all vectors are trivially eigenvectors. 

have solutions as long as the discriminant 


For /3 > 7 , if / or <7 is odd, eqs. 
of M, 5' = — eh + gf, admits square roots modulo 2^“"^; we are led to an 

analysis similar to the one above concerning D = 2^'^P. 

If both / and g are even, then S = (±1 -|- 2"<)l2 mod 2'^~^^. This case implies 
a = 27 , so (5 = 5'. All vectors modulo are eigenvectors. For /3 > 7 -|- 2, eq.(^) 
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can be divided by a factor 2 , leading to: 

Fk^ + Ek — G = 0 mod or the analogous equation in k', 

where F = f/2, E = {e — h)/2, G = gj2. (91) 

By definition, M 7 ^ 0 mod 2, so the equation ( ^ 8 |) implies that e, h, and E are 
odd. Therefore, since E is the derivative of the polynomial ( pl| ) modulo 2, lemma 
|I] associates to each solution k modulo 2 a unique solution modulo for any 

(3. The analysis modulo 2 is straightforward. 


(F,G) = (0,0) 
(F,G) = (1,0) 
(F,G) = (0,1) 
(F,G) = (1,1) 


A; = 0 


or k! = ^ 


k = 0 or A: = 1 
k' = 0 or A:^ = 1 


no solution 


(92) 


Equivalently, in this last case (/, g even) there are 2'’'^^ solutions for any large (3 if 
5 = 1 mod 8 (the first three equations), and none if 5 = 5 mod 8 (the last equation), 
which are the only possible values of 5 in this case. 
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Figure captions 


Fig. 1. Husimi function of the basis state |g 2 )Af=ii,K= 0 ! plotted both in linear 
scale (gray levels; the high values are dark) and logarithmic (contour) scale. The 
linear scale shows a concentration of the density along the line {q = q 2 = 2 / 11 }, 
whereas the logarithmic plot shows the vertical string of zeros along the line {q = 
q 2 + 1 / 2 }. The function has the period xjWypl in the variable z = {q — ip)/V2. 


Fig. 2. For the three admissible values ofA^ = ll, A^ = 33 = llx3, N = 


121 = 11 ^, we plot free eigenmodules An of S'fig = 


generating vectors are respectively vn = 


fig 

1 

6 


2 1 
3 2 


I, rescaled by 1/A^. The 


, ^'33 = 


1 

6 


■^121 = 


1 

94 


We notice that the string An/II is a sublattice of A 33/33 (cf. eq. (^)), and of 
A 121 /I 2 I (cf. lemma H). 


Fig. 3. Husimi functions of two k = 0 pure theta eigenstates of for 

N = 33 (left) and N = 121 (right); we show the linear density scale using gray 
levels and indicate the position of the zeros by tiny circles. In both cases, the zero 
patterns are the rescaled affine strings A'^/Ny/2, where A)y are translates of the 
eigenmodules A 33 , resp. A 121 of figure 2. For both eigenstates, the eigenangle is 
— 7 r/ 4 , of respective degeneracies 2 and 3. 


Fig. 4. Husimi function of the k = 0 eigenstate of ClSfjg associated with the 
invariant square sublattice A(q, for = 121 (same representation as in fig. 3). The 
eigenangle is — vr/4, hence this state is degenerate with the pure theta eigenstate 
shown in fig. 3 (right). 


Fig. 5. Husimi functions of two k = 0 cristalline eigenstates of Hsflg, for N = 33. 
In both cases, we use translates of the eigenmodule A 33 of figure 2 (center). On the 
left, the coset Oi contains li = 11, its period is pi = 2; the eigenangle is 2tt/3 — 7 r /4 
(triply degenerate), and the crystal has the period (—1 + 6i)/ll\/2. On the right, 
O 2 contains I 2 = 3, its period is p 2 = 10; the eigenangle is —27r/10 — 7 r /4 (non¬ 
degenerate), and the crystal has the period l/3-v/2. The larger circles materialize 
double zeros. 


Fig. 6 . Husimi functions of two k = 0 cristalline eigenstates of Hg^g. On the left, 
for the square value = 81, we use a cycle of translates of the lattice A(,q, starting 
from = (3 -|- 3i)/9; this cycle is of period 6 , the eigenangle is 2it/6 — 7r/36 (non¬ 
degenerate), and the zero pattern inherits a square periodicity of side l/3\/2- On 
the right, the value A^ = 55 is not admissible, we use a cycle of rotations starting 
from a linear string A” containing the invariant sublattice 5Aii (An is shown in 
figure 2 (left)). The cycle is of period prot = 3, yielding an eigenangle 2tt/3 — 3'k/A 
(four times degenerate), and the crystal has the period (—1 -|- 6i)/ll\/2. Moreover, 
the function vanishes on the string An/ll\/ 2 . 

Fig. 7. Husimi function of two k = 0 eigenstates of Usf^^, for N = 33 (left) and 
N = 59 (right), both being admissible values. These state are built as described 
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in subsection (O^ i-e- using translates of invariant strings. On the left, we use a 
translation cycle O 3 of period 10 containing ^3 = 1, the eigenangle is 27r/30 — vr/d 
(doubly degenerate); on the right, the cycle is of period 58, and the eigenangle 
7 r /4 — 27r/29 is non-degenerate. In both cases, the state presents no particular 
periodicity, but rather a ‘scar’ at a classical fixed point ((1 -|- i)/ 2\/2 on the left, 0 
on the right). Nevertheless, the zero patterns seem well-distributed on the torus. 
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